Keywords 

Bottleneck, conservation genetics, 
demographic survey, Parachondrostoma 
toxostoma, rivers, species distribution 
models, temporal trends. 

Correspondence 

Ivan Paz-Vinas, Laboratoire Evolution et 
Diversite Biologique (EDB), UMR 5174 
(CNRS - UPS - ENFA), 118 route de 
Narbonne, F-31062 Toulouse cedex 4, 
France. Tel: (+33) 5 61 55 67 47; Fax: (+33) 
5 61 55 73 27; E-mail: ivanpaz23@gmail.com 

Funding Information 

This study is part of the European project 
"IMPACT". This project has been carried out 
with financial support from the Commission 
of the European Communities, specific RTD 
programme "IWRMNET." 



Ecology and Evolution 



Combining genetic and demographic data for prioritizing 
conservation actions: insights from a threatened fish 
species 

Ivan Paz-Vinas 1 ' 2 ' 3 , Lise Comte 1,2 , Mathieu Chevalier 1,2,4 , Vincent Dubut 5 , Charlotte Veyssiere 1,2 , 
Gael Grenouillet 1,2 , Geraldine Loot 3,2 & Simon Blanchet 3,1 

1 UMR5174 EDB (Laboratoire Evolution & Diversite Biologique), Centre National de la Recherche Scientifique (CNRS), Ecole Nationale de Formation 
Agronomique (ENFA), Universite Paul Sabatier, 118 route de Narbonne, F-31062, Toulouse Cedex 4, France 
2 UMR 5174 (EDB), UPS, Universite de Toulouse, 1 18 route de Narbonne, F-31062, Toulouse Cedex, France 

3 Centre National de la Recherche Scientifique (CNRS), Station d'Ecologie Experimental du CNRS a Moulis, USR 2936, F-09200, Moulis, France 
4 UMR 5245 EcoLab (Laboratoire Ecologie Fonctionnelle et Environnement), CNRS, F-31062, Toulouse, France 

5 IMBE - UMR 7263, Aix-Marseille Universite, CNRS, IRD, Centre Saint-Charles, Case 36, 3 place Victor Hugo, F-13331, Marseille Cedex 3, France 



Received: 17 January 2013; Revised: 16 May 



Abstract 

Prioritizing and making efficient conservation plans for threatened populations 
requires information at both evolutionary and ecological timescales. Neverthe- 
less, few studies integrate multidisciplinary approaches, mainly because of the 
difficulty for conservationists to assess simultaneously the evolutionary and eco- 
logical status of populations. Here, we sought to demonstrate how combining 
genetic and demographic analyses allows prioritizing and initiating conservation 
plans. To do so, we combined snapshot microsatellite data and a 30-year-long 
demographic survey on a threatened freshwater fish species {Parachondrostoma 
toxostoma) at the river basin scale. Our results revealed low levels of genetic 
diversity and weak effective population sizes (<63 individuals) in all popula- 
tions. We further detected severe bottlenecks dating back to the last centuries 
(200-800 years ago), which may explain the differentiation of certain popula- 
tions. The demographic survey revealed a general decrease in the spatial 
distribution and abundance of P. toxostoma over the last three decades. We 
conclude that demo-genetic approaches are essential for (1) identifying popula- 
tions for which both evolutionary and ecological extinction risks are high; and 
(2) proposing conservation plans targeted toward these at risk populations, 
and accounting for the evolutionary history of populations. We suggest that 
demo-genetic approaches should be the norm in conservation practices. 
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Introduction 

Prioritizing and making appropriate plans to manage 
and conserve threatened species is a complex task. 
Global changes simultaneously affect multiple facets of 
individual species, making predictions difficult 
(Margules and Pressey 2000; McMahon et al. 2011). For 
instance, global changes such as habitat fragmentation 
or climate change can affect the genetic diversity (Olivieri 



et al. 2008; Blanchet et al. 2010), the demographic 
dynamics (Julliard et al. 2004; Dunham et al. 2008), the 
evolution of life-history traits (Conover et al. 2009; 
Blanchet and Dubut 2012), and/or the spatial distribu- 
tion of species (Parmesan 2006; Buisson et al. 2008). 
Accordingly, the conservation biologists' toolbox 
includes several methods which emerged from multiple 
disciplines such as population genetics, population ecol- 
ogy, and biostatistics (Guisan and Zimmermann 2000; 
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Green et al. 2005; Excoffier and Heckel 2006). Neverthe- 
less, most conservation studies focus on a single facet of 
species health (e.g., the genetic diversity), and hence 
provide only partial information for biodiversity man- 
agement and conservation (Frankham 2010; Geist 2011; 
Loss et al. 2011). 

Integrative studies are, however, increasingly acknowl- 
edged as being valuable from a conservation standpoint 
(Purvis and Hector 2000; Geist 2011; Loss et al. 2011). 
For instance, at the community level, Devictor et al. 
(2010) showed that there was a strong spatial mismatch 
between phylogenetic, functional, and taxonomic mea- 
sures of bird biodiversity. These measures provide differ- 
ent but complementary information, suggesting that 
reserve designs should be optimized accordingly (Devictor 
et al. 2010). Similarly, at the population level, diverse 
measures classically used to assess the health of a popula- 
tion (e.g., effective population size, abundance, and dis- 
persal rate) provide complementary information that 
should be integrated into common analyses to set efficient 
conservation plans (e.g., Osborne et al. 2010, 2012). For 
instance, demographic monitoring programs (hereafter, 
DMPs) provide useful information regarding the ecologi- 
cal status of populations and enable predictions on future 
distributions under global change scenarios, whereas pop- 
ulation genetics studies (hereafter, PGSs) obtain informa- 
tion regarding the evolutionary status of populations and 
their potential resistance to rapid environmental changes 
(Smith and Bernatchez 2008). Because evolutionary and 
ecological timescales and processes are sometimes con- 
founded (Carroll et al. 2007), it is of prime importance 
to merge evolutionary and ecological information to (1) 
identify the populations that need to be prioritized for 
conservation actions; and (2) implement effective 
long-term management and conservation of endangered 
populations (Osborne et al. 2012). 

The use of population genetics in biodiversity conser- 
vation has increased considerably in the last decades 
(Frankham 2010). Low genetic diversity in natural popu- 
lations has been generally associated with pervasive 
effects such as inbreeding depression, loss of evolution- 
ary potential, and the accumulation of deleterious muta- 
tions (Saccheri et al. 1998; Frankham 2010). These 
effects theoretically increase extinction risks, and are 
expected to be stronger in populations under anthropo- 
genic or natural stresses (Spielman 2004). Accordingly, 
PGSs generally aim at (1) describing the genetic status 
of populations (i.e., genetic diversity and structure 
assessed during a snapshot survey, Schwartz et al. 2007); 
(2) identifying historical and contemporary factors 
affecting the genetic diversity of populations (Manel 
et al. 2003; Dubut et al. 2012); and (3) inferring past 
and contemporary demographic parameters such as 



effective population sizes (N e ) (Storz and Beaumont 
2002). Although PGSs provide key information about 
demographic processes, linking genetics and population 
demography remains tricky (Osborne et al. 2012). For 
instance, the link between N e and census population size 
(N c ) is notoriously difficult to assess (Luikart et al. 
2010; Belmar-Lucero et al. 2012; Palstra and Fraser 
2012), and genetic bottlenecks (i.e., strong decreases in 
AT e ) can be detected even in the absence of demographic 
bottlenecks (Broquet et al. 2010; Chikhi et al. 2010). 
Furthermore, the effects of particular threats may be 
undetected through PGSs due to the lag time that often 
exists between an ecological cause and its evolutionary 
consequence (Landguth et al. 2010). 

Analyses based on demographic data can overcome 
some of these gaps (Nichols and Williams 2006; Linden- 
mayer et al. 2010). DMPs provide information about the 
current status of populations by allowing the inference 
of key demographic parameters such as abundance and/ 
or occurrence (Royle and Dorazio 2006). Combined 
with time series analyses, DMPs also permit the investi- 
gation of temporal trends and hence the identification 
of the causes and consequences of population declines 
or changes in spatial distribution (Daufresne et al. 
2004). Additionally, these surveys are useful for the early 
detection of the effects of threats on populations as well 
as "ecological surprises" (Doak et al. 2008), which is 
notoriously difficult using only PGSs (Julliard et al. 
2004; Lindenmayer et al. 2010). Finally, long-term and 
large spatial-scale surveys are of prime interest and may 
allow predictions about the future status of populations 
in a changing world through the use of species distribu- 
tion models for instance (Guisan and Zimmermann 
2000). 

In this study, we attempt to demonstrate how combin- 
ing PGSs and DMPs provides baseline information for 
prioritizing and initiating management and conservation 
plans. We focused on an endangered freshwater fish spe- 
cies (i.e., the South-west European nase Parachondrostoma 
toxostoma, Vallot 1837) which is considered vulnerable 
throughout its restricted native range (i.e., Southern 
France, Crivelli 2006). We used a microsatellite dataset 
gathered at the river basin scale (i.e., the Garonne river 
basin, South- Western France) to (1) describe the genetic 
diversity and structure of P. toxostoma populations, and 
(2) detect and quantify both contemporary and past N e 
(i.e., contraction or reduction in N e over time), as well as 
to date main changes in N e following the last glacial max- 
imum (i.e., approximately 10,000 years ago). In parallel, 
we used a demographic survey performed at the same 
spatial scale over the last three decades to (3) identify 
temporal trends in species abundance at the Garonne 
river basin scale; and (4) assess the current spatial distri- 
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bution of the species and changes in the distribution over 
the last three decades. 

Materials and Methods 
Biological model 

Parachondrostoma toxostoma is a threatened freshwater 
fish species of the Cyprinidae family endemic to France 
and Switzerland, where its native range area is restricted 
to the Rhone, Adour and Garonne river basins. This 
species is listed as vulnerable in the IUCN red list, in 
the Annex II of the European Union Habitats Directive 
and in Appendix III of the Bern Convention (Crivelli 
2006). The range of the species has been strongly 
reduced due to water pollution, habitat fragmentation 
by dams and weirs, artificial water releases and hybrid- 
ization with a nonnative species, Chondrostoma nasus 
(Costedoat et al. 2007). Our study focuses on the Gar- 
onne river basin, which hosts the major stock of pure 
P. toxostoma (i.e., not introgressed by the C. nasus gen- 
ome). This highlights the urge for conservation actions 
directed toward the Garonne drainage in order to pre- 
serve the P. toxostoma species. 



Population genetics study 

Sampling design 

Ninety-two sampling sites belonging to 34 rivers of the 
Garonne river basin were investigated using electrofishing 
in 2010 and 2011 (Fig. SI). We did not catch P. toxos- 
toma at 76 sites. Two hundred and 30 individuals of 
P. toxostoma were sampled at sixteen sites (Table 1, 
Fig. 1). Thus, we assume that these sixteen sites are 
representative of the current P. toxostoma populations. 
However, due to the low numbers of individuals captured 
at some sampling sites, individuals from sites belonging 
to the same river were pooled for subsequent analyses. All 
genetic analyses were therefore conducted at the river 
level (mriver = 9). A small fragment of pelvic fin was col- 
lected and stored in 90% ethanol. Individuals were all 
released alive at their sampling site. 

Genotyping 

We used a salt-extraction protocol to extract genomic 
DNA from pelvic fins (Aljanabi and Martinez 1997). 
Fifteen microsatellite loci previously developed and/or 



Table 1. Parachondrostoma toxostoma sampling sites information. 



River 


Code 


Location 


Latitude 


Longitude 


PGS 


W(PGS) 


DMP 


V(DMP) 


ARRATS 


ARR 


Aubiet 


N 43°38'48" 


E 0°46'4S" 






X 


13 


AUROUE 


AUR 


L'isle-Bouzon 


N 43°54'32" 


E 0°43'45" 






X 


13 


AVEYRON 


AVE 


Feneyrols 


N 44°07'52" 


E 1°48'51" 


X 


5 










Monteils 


N 44° 17'09" 


E 2°00'07" 


X 


4 






ARIEGE 


ARI 


Venerque 


N 43°26'13" 


E 1-26'15" 






X 


8 


PETITE BARGUELONNE 


BAR 


Montbarla 


N 44° 12'34" 


E 1°03'40" 


X 


9 


X 


17 


CELE 


CEL 


Boussac 


N 44°35'46" 


E 1°55'02" 


X 


7 










Sainte Eulalie 


N 44°35'36" 


E 1°52'25" 


X 


8 










Sauliac-sur-Cele 


N 44°31 '09" 


E 1°42'58" 


X 


25 


X 


1 1 


COUZE 


COU 


Bayac 


N 44°48'1 6" 


E 0°43'4S" 






X 


14 


ELLE 


ELL 


Terrason-Lavilledieu 


N 4S°08'51" 


E 1°15'37" 


X 


25 






GARONNE 


GAR 


Muret 


N 43°27'36" 


E 1°19'52" 






X 


10 


HERS 


HER 


Besset 


N 43°05'03" 


E 1°50'24" 


X 


4 


X 


10 






Calmont 


N 43°17'10" 


E 1°37'S9" 


X 


25 






LOUGE 


LOU 


Fousseret 


N 43°16'27" 


E 1"04'07" 


X 


8 


X 


13 


SALAT 


SAL 


Touille 


N 43°04'38" 


E 0°B8'0S" 


X 


25 






SAVE 


SAV 


Espaon 


N 43°25'20" 


E 0°51'21" 


X 


18 






VENDINELLE 


VEN 


La Salvetat Lauragais 


N 43°32'22" 


E 1°48'15" 






X 


18 


VERE 


VER 


Cahuzac-sur-Vere 


N 43°59'12" 


E 1°53'43" 






X 


17 


VIAUR 


VIA 


La Calquiere 


N 44°09'12" 


E 2°12'1S" 


X 


13 










Saint Just 


N 44°07'24" 


E 2°21'S7" 


X 


23 










Navech 


N 44°09'25" 


E 2°23'18" 


X 


25 










Serres 


N 44°12'29" 


E 2°31'25" 


X 


6 






VOLP 


VOL 


Plan 


N 43°10'16" 


E 1°07'07" 






X 


8 



PGS (for Point Genetic Study) indicates whether the site has (X) or not (-) been sampled for genetic analyses. W<pgs) indicates the number of indi- 
viduals sampled per site for genetic analyses. DMP (for Demographic Monitoring Program) indicates whether the site has (X) or not (-) been 
selected for analyses of temporal trends in abundance. V( DMP ) indicates the number of years considered in the time series. 
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Figure 1. Map of the Garonne river basin (South-Western France) 
representing (1) sites where Parachondrostoma toxostoma was 
sampled for the genetic analyses (green circles) and (2) sites that have 
been selected for analyses of temporal trends in population 
abundances (black triangles). 

evaluated for P. toxostoma (Dubut et al. 2010) were coam- 
plified using two multiplexed polymerase chain reactions 
(PCRs; see Table SI for details on loci and primers concen- 
trations). PCR amplifications were performed with 
5-20 ng of genomic DNA and using the QIAGEN® Multi- 
plex PCR Kit (Qiagen, Valencia, CA). PCRs were carried 
out under conditions described by Dubut et al. (2010). 
Genotyping was performed on an ABI PRISM™ 3730 
Automated Capillary Sequencer (Applied Biosystems, 
Foster City, CA) at the "Genopole Toulouse Midi- 
Pyrenees" (France). Allele sizes were scored using the soft- 
ware GENEMAPPER® v.4.0 (Applied Biosystems). 

Descriptive genetic analyses 

The presence/ absence of large allele dropouts and null 
alleles was determined using the software MICRO - 
CHECKER 2.3 (Van Oosterhout et al. 2004). Departures 
from Hardy-Weinberg (HW) equilibrium were estimated 
using the program GENEPOP v4.0 (Rousset 2008). Levels 
of significance for HW were adjusted using the false dis- 
covery rate (FDR) procedure (Benjamini and Hochberg 
1995). Linkage disequilibrium among loci within sites was 
tested with the program FSTAT 2.9.3.2 (Goudet 1995). 

The mean number of alleles per site, the average 
observed (H 0 ) and expected (H e ) heterozygosity over loci, 
as well as H Q and H e per loci per site were estimated 
using ARLEQUIN 3.5 (Excoffier and Lischer 2010). 
We used a rarefaction procedure, as implemented in the 



software ADZE 1.0 (Szpiech et al. 2008), to estimate 
allelic richness (Petit et al. 1998) for each site, considering 
minimum sample sizes of N = 8 and N = 18 individuals. 

Population structure 

A Bayesian model-based clustering approach was used to 
search for the occurrence of independent genetic groups 
(i.e., clusters, K) in our dataset (as implemented in 
STRUCTURE 2.3.3; Pritchard et al. 2000; Falush et al. 
2003, 2007; Hubisz et al. 2009). The burn-in length of the 
Markov Chain Monte Carlo (MCMC) was set to 50,000 
followed by 200,000 iterations. The admixture model and 
the correlated allele frequencies model were used with 
priors on population sampling location (Hubisz et al. 
2009). Ten runs were conducted for each K value, with K 
ranging from 1 to 10. We used CORRSIEVE 1.6.2 (Cam- 
pana et al. 2011) to combine two approaches aiming at 
determining K: the AK test (Evanno et al. 2005) and the 
AF st test (Campana et al. 2011). 

To further assess the levels of genetic differentiation 
among P. toxostoma sites, two different indices were esti- 
mated: pairwise F st (Weir and Hill 2002) and the unbiased 
pairwise D est (Jost 2008), calculated using ARLEQUIN 3.5 
and SMOGD (Crawford 2010), respectively. 

Demographic history inference and current 
N e estimation 

We used two different approaches for inferring past 
changes in the effective population size (i.e., expansions 
or contractions) of P. toxostoma. 

The first method, implemented in the BOTTLENECK 
vl.2.02 software (Cornuet and Luikart 1996; Piry 1999), 
uses summary statistics of the genetic diversity to assess sig- 
nificant deviations from mutation/drift equilibrium. Signifi- 
cant heterozygosity excesses are considered as evidence of 
recent bottlenecks, whereas significant heterozygosity defi- 
ciencies can be interpreted as signals of recent population 
expansion (Luikart and Cornuet 1998). We performed 
analyses considering two different microsatellite evolution 
models: the stepwise mutation model (SMM) and the two- 
phase model (TPM). For the latter, we set the percentage of 
multistep mutations at 30%. We tested the significance of 
mutation/drift equilibrium deviations for the two models 
using Wilcoxon's signed rank tests. To account for multiple 
comparisons, we applied the FDR procedure (Benjamini 
and Hochberg 1995). The second method is the full-likeli- 
hood Bayesian approach implemented in the program 
MSVAR 1.3 (Beaumont 1999; Storz and Beaumont 2002). 
This coalescent-based method relies on a hierarchical Bayes- 
ian model to detect, date, and quantify past demographic 
changes. The model assumes that a stable, closed population 
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of ancestral size Nj increased or decreased exponentially to 
its current size N 0 (i.e., its current N e ) over a time interval 
of T a years. This method uses all the information contained 
in the data and lognormal priors to infer the parameters of 
the model <I> = {N 0 , N u T a , 9}, where 9 = 4N 0 fi and ft is 
the mutation rate. The posterior probability density of <J> is 
assessed via MCMC algorithms. Microsatellite loci are 
assumed to be independent and to evolve under a strict 
SMM. For each river-scale analysis, we performed four inde- 
pendent runs of 5 x 10 9 steps, considering different starting 
values and means for priors and hyperpriors for each run 
(Goossens et al. 2006). We set a generation time of 3 years 
for P. toxostoma (Keith et al. 2011). Parameters were 
thinned with an interval of 5 x 10 4 steps, resulting in out- 
put files with 1 x 10 5 values. We discarded the first 10% of 
the chains as burn-in to prevent bias induced by the starting 
values on parameter estimation. The convergence of the 
MCMC chains was checked with the Gelman and Rubin 
analysis implemented in the R package CODA (Gelman and 
Rubin 1992; Plummer et al. 2006). For each analysis, pos- 
terior parameter values obtained by the four independent 
runs were pooled together and subsequently used to calcu- 
late the median and the 5-95% quartiles for N 0 , Nj, and T a . 
We also calculated these statistics for the ratio log 10 (N 0 /N 1 ). 
Negative values of this ratio indicate that the population has 
experienced a decrease in effective population size, while 
positive values characterize demographic expansions. This 
approach was also used to estimate a current N e at the Gar- 
onne river basin scale. To do so, we ran MSVAR by pooling 
all individuals from all rivers in a single analysis. At such a 
scale, estimates of current N e were compared to those esti- 
mated using the linkage disequilibrium-based approach 
implemented in LDNe (Waples and Do 2008). LDNe was 
not used at the river scale due to its propensity to give nega- 
tive N e estimates (which are interpreted as infinity estimates, 
Waples and Do 2008) for most rivers. MSVAR 1.3 runs 
were performed on an ALTIX ICE 8200 EX and UV com- 
puter cluster (Silicon Graphics International, Fremont, CA) 
hosted by the CALMIP group at the University Paul Saba- 
tier (Toulouse, France). 

Demographic monitoring data 

Database description 

We used the surveillance monitoring database of the 
French National Agency for Water and Aquatic Environ- 
ments (i.e., ONEMA) to carry out demographic trend 
and species distribution analyses. This database includes 
an extensive spatiotemporal set of monitoring surveys of 
French freshwater fish populations, representative of all 
fish assemblages and covering varying degrees of anthro- 
pogenic disturbances (Poulet et al. 2011). Surveys were 



conducted according to standard electrofishing procedures 
(Poulet et al. 2011). We used this database to (1) identify 
temporal trends in population abundance of P. toxostoma 
at 12 sampling locations; (2) assess the current spatial dis- 
tribution of this species in the Garonne river basin; and 
(3) investigate whether the spatial distribution of this spe- 
cies in the Garonne river basin has declined or expanded 
over the last three decades. 

Temporal trends in abundance 

From this dataset, we selected all sites belonging to the 
Garonne river basin that have been sampled and investi- 
gated for P. toxostoma abundance for at least 8 years. This 
resulted in the selection of twelve sites (Table 1, Fig. 1) 
for which time series ranged between 8 and 18 years and 
occurred between 1991 and 2010. As sampling procedures 
were standardized over years, abundances (expressed as 
the number of individuals per m 2 ) were directly compara- 
ble across years. It is noteworthy that (1) this database 
and the genetic database have been gathered during inde- 
pendent research projects; and (2) P. toxostoma is rela- 
tively rare in this area (Fig. SI), which both explain why 
demographic and genetic data are not available for all 
sites (see Table 1). Some sites for which long-term demo- 
graphic data were available have been unsuccessfully sam- 
pled for genetic, and inversely, some sites where genetic 
data were available had time series that were not long 
enough to be analyzed (i.e., <8 years). 

First, we assessed the strength and significance of tempo- 
ral trends at these sites, by using a modified Mann-Kendall 
trend test that we independently applied to each time series 
(Hamed and Rao 1998). In this test, the Mann-Kendall's S 
statistic (Kendall 1962) provide an estimate of the strength 
of the association between time and the response variable, 
while accounting for temporal autocorrelation present in a 
time series (Hamed and Rao 1998). 

Second, we assessed whether or not these twelve time 
series showed an overall significant trend. For this pur- 
pose, we performed a meta-analysis (Gurevitch and 
Hedges 1993) on the twelve Mann-Kendall's trend statis- 
tics S calculated in the first step. We applied a mixed lin- 
ear model approach using maximum likelihood, in which 
we assumed that the 12 time series included in the 
meta-analysis share a common effect size with a random 
variation among the twelve time series. 

Current spatial distribution and recent 
distribution changes 

We used the database described above to assess changes 
in the spatial distribution of P. toxostoma on the Garonne 
river basin over two distinct periods, separated by a time 
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span of 10 years (i.e., "past period": 1980-1992, and "cur- 
rent period": 2003-2009). To account for potential sam- 
pling bias when comparing spatial distributions over time 
based on datasets not originally collected for this purpose 
(Shaffer et al. 1998; Shoo et al. 2006), we modeled the 
spatial distribution of the species across the French 
hydrographic network as a function of several climatic 
and environmental variables. 

Accurately modeling species distribution requires per- 
forming analyses at the entire species range scale, so as to 
encompass all environmental conditions (Austin 2007). 
Therefore, for both time periods, initial models were cali- 
brated at the French scale. We selected 3549 sites sampled 
over the 1980-1992 period and 3543 sites sampled over 
the 2003-2009 period scattered across France (see Fig. 
S2). The occurrence of the species was modeled indepen- 
dently for both time periods as a function of habitat and 
climatic data strongly related to fish spatial distributions 
(Buisson et al. 2008): elevation (m), slope (%), 
upstream-downstream position (G), mean temperature of 
the coldest quarter (°C), mean temperature of the warm- 
est quarter (°C), temperature variability, cumulated 
precipitations of the wettest quarter (mm), cumulated 
precipitations of the driest quarter (mm), and precipita- 
tion variability (Hijmans et al. 2005). 

To account for uncertainty in estimating species range, 
we used a modeling approach allowing us to produce 
maps of species habitat suitability (e.g., Puschendorf et al. 
2009; Grenouillet et al. 2011). Specifically, we used an 
ensemble modeling approach based on a consensus model 
averaging the probabilities of occurrence predicted by 
eight single-species distribution models (Marmion et al. 
2009), as well as three threshold setting methods allowing 
the conversion of occurrence probabilities into binary 
data (i.e., presence or absence, Liu et al. 2005), and 30 
iterations (see Appendix SI for details on models' 
implementation) . 

The calibrated models set at the French scale were then 
used to predict the binary predictions of occurrence of 
the species for the two distinct periods in the hydro- 
graphic network of the Garonne river basin. The spatial 
distribution of the species for each time period was calcu- 
lated as the length of the hydrographic network occupied 
by the species (e.g., Fagan 2002) in the Garonne river 
basin (expressed in % of the total network length). How- 
ever, because the ability to detect changes in the spatial 
distribution of species may be confounded by the uncer- 
tainty arising from methodological strategies (e.g., thresh- 
old effect, Nenzen and Araujo 2011), temporal changes in 
the occupied stream length were evaluated using a linear 
model that controlled for the threshold effect. A linear 
model was thus fitted to the spatial distribution of P. tox- 
ostoma in both periods where the threshold-setting 



method and the period were used as explanatory vari- 
ables. The change (i.e., extension or contraction) was then 
provided by the least-squares means intercepts of the con- 
temporary period-group effect. Temporal trends analyses 
and spatial distribution models have been developed 
under the R statistical software 2.13.0 (R Development 
Core Team 2011). 

Results 

Population genetics study 

Descriptive genetic analyses 

After applying the FDR controlling procedure, no null 
alleles were detected in our dataset, there were no signifi- 
cant deviations from HW for any loci or any population 
(Tables S2 and S3), and we failed to detect significant 
linkage disequilibrium between pairs of loci (Table S4). 

Overall, genetic diversity estimates were low (Fig. 2A 
and B, Table S3). Loci were weakly polymorphic at the 
basin scale (2-6 alleles per locus), with some loci being 
monomorphic at the river scale (na = 1; Table S2). Aver- 
age H e and H G values across loci within rivers were mod- 
erately low (H e = 0.320-0.450; H G = 0.315-0.482), as well 
as mean number of alleles and allelic richness estimates 
(AR 8 = 1.868-2.536 alleles per river; AR 18 = 2.147-3.037; 
Fig. 2A and B, Table S3). It is noteworthy that the 
Save River (SAV) displayed the lowest genetic diversity 
estimates (Fig. 2A and B, Table S3). 

Population structure 

The ten runs of the Bayesian clustering analysis were con- 
vergent. The AK and AF st tests revealed three distinct 
clusters K = 3 (Fig. 3A-B). Most of the populations were 
hardly differentiable and were characterized by the occur- 
rence of a main cluster, whose frequency range was from 
62% (CEL) to 98% (VIA). Only SAV and HER were dis- 
criminated from the rest of the Garonne river basin, each 
site corresponding to a distinct cluster (Fig. 3C). Overall, 
genetic differentiation values between rivers were weak to 
moderate and ranged between 0.003 and 0.244 and 0.003 
and 0.281 for F st and D estj respectively (Table 2). All but 
five pairwise F st values were significant (Table 2). The 
stronger differentiations were found between SAV/VIA 
(F s , = 0.244; D est = 0.097) and SAV/BAR (F st = 0.117; 
D est = 0.281). 

Demographic history inference and current N e 
estimation 

According to the BOTTLENECK software, and after 
corrections for multiple tests, there was no significant 
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Figure 2. Maps representing (A) the allelic richness per population considering a minimum sample size of 8 (color scale), (B) the expected 
heterozygosity per population (color scale), (C) the past effective population size (N-i; left number in the bubbles, see also Table S6), the current 
effective population size (W 0 ; right number in the bubbles, see also Table S6), the time of the beginning of the bottlenecks (in years backward in 
time; numbers in brackets, see also Table S6), and the magnitude of bottlenecks (i.e., Log 10 (W0/W1): color scale, see also Table 56), and (D) the 
value of the Mann-Kendall's S statistic (color scale) and the significance of Mann-Kendall trend tests for each time series: Asterisks (*) denote 
significant (i.e., P < 0.05) temporal trends. For all panels, the three-letter code in each bubble corresponds to the river codes (see Table 1). 



evidence for demographic changes in the Garonne river 
basin (Table S5). On the contrary, the MSVAR analyses 
revealed significant signals of bottleneck in all rivers 
(Fig. 2C, Table S6). The magnitude of these bottlenecks, 
as indicated by the median values of the log 10 {N 0 /Ni) 
ratio, ranged between -0.705 (ELL) and -1.345 (HER; 
Fig. 2C, Table S6). Overall, N 0 estimates (i.e., the current 
N e of populations) were similar across rivers, with 



medians ranging from 7 (HER) to 63 individuals (SAL). 
Concerning ancestral population sizes (N^), median val- 
ues ranged from 5286 (LOU) to 9155 individuals (HER; 
Fig. 2C, Table S6). These bottlenecks were estimated to 
have occurred between 192 (HER) and 727 years ago 
(SAL). The MSVAR method has often been considered as 
more powerful than the BOTTLENECK method 
(Williamson-Natesan 2005; Girod et al. 2011), which 
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ELL CEL £ ^ VIA g SAL HER SAV 

CD < _i 

Figure 3. Analysis of the population structure of Parachondrostoma toxostoma in the Garonne river basin. (A) and (B) represent the results from 
AK and AF st tests, respectively. (C) is a barplot representing the results of the Bayesian clustering analysis of microsatellites using STRUCTURE for 
K = 3. 



Table 2. Population pairwise F sl (upper half-matrix) and pairwise D est (lower half-matrix) values calculated between all rivers (denoted by their 



three-level code). 


Code 


AVE 


BAR 


CEL 


ELL 


HER 


LOU 


SAL 


SAV 


VIA 


AVE 




0.117 


0.067 


0.070 


0.042 


0.013 ns 


0.014 


0.109 


0.005 ns 


BAR 


0.056 




0.102 


0.052 


0.054 


0.025 


0.026 


0.130 


0.017 


CEL 


0.018 


0.031 




0.023 


0.023 


0.003 ns 


0.012 


0.089 


0.010 


ELL 


0.035 


0.008 


0.003 




0.032 


0.008 ns 


0.014 


0.115 


0.008 


HER 


0.132 


0.165 


0.077 


0.096 




0.069 


0.068 


0.077 


0.122 


LOU 


0.054 


0.057 


0.029 


0.029 


0.019 




0.004 ns 


0.114 


0.050 


SAL 


0.049 


0.096 


0.033 


0.034 


0.013 


0.024 




0.090 


0.037 


SAV 


0.262 


0.281 


0.221 


0.230 


0.241 


0.265 


0.228 




0.244 


VIA 


0.026 


0.093 


0.044 


0.034 


0.033 


0.015 


0.010 


0.097 





For pairwise F st , significant values at level 0.05 after false discovery rate (FDR) correction are in bold. Nonsignificant pairwise F st are denoted by 
ns. 
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may explain the discrepancy observed between these two 
methods. 

The analysis performed at the Garonne river scale con- 
firmed the low estimates of current N e found at the river 
scale. Indeed, at this scale, MSVAR provided an estimate 
of 147 individuals (5-95% quartiles: 35.6-534.4) in the 
whole drainage, whereas LDNe provided a global estimate 
of 74.6 individuals (95% CI: 54.4-104.6). 

Demographic monitoring data 

Temporal trends in abundance 

Five out of the twelve populations (i.e., HER, VEN, AUR, 
CEL, VER) showed a significant negative trend (P < 0.05; 
S < 0), one population (COU) showed a significant posi- 
tive trend (P < 0.01; S = 23) whereas the remaining six 
populations (VOL, LOU, ARI, GAR, ARR, BAR) showed 
no significant trend in abundance (Fig. 2D, Table S7). 
Overall, the mixed model meta-analysis revealed a signifi- 
cant (P < 0.001) negative trend indicating a global 
decrease in the abundance of P. toxostoma populations in 
the Garonne river basin. 

Modeling species distribution 

The stream length occupied by the species was estimated 
at 24.0% (±2.5 SE) of the total river basin stream length 
in 1980-1992 (Fig. 4A) and 20.9% (±2.6 SE) in 2003- 
2009 (Fig. 4B). This represented an overall decrease of 
3.2% (P < 0.01) with respect to the whole river basin, 
and of 13.1% of P. toxostoma's 1980-1992 distribution 
(Figs. 4C, 5). The habitat suitability for the species 
decreased in the middle part of the river basin between 
1980 and 1992 and 2003 and 2009 periods (Fig. 4). 

Discussion 

What did we learn from genetic data? 

Using a full-likelihood Bayesian approach (as imple- 
mented in MSVAR, Storz and Beaumont 2002), we 
showed that all P. toxostoma populations have experi- 
enced significant decreases in effective population size 
(N e ), with reductions of more than 99% of their prebot- 
tleneck long-term N e . We further showed that: (1) in all 
populations, bottlenecks started 192-727 years ago, and 
are hence relatively recent (i.e., within the last millen- 
nium); and (2) all populations show extremely low 
current N e . Attempting to identify the causes of such bot- 
tlenecks would be highly speculative without further data 
and analyses. If natural causes (climatic or hydrological 
shifts) cannot be ruled out, anthropogenic causes are also 
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Figure 4. Spatial distributions of Parachondrostoma toxostoma 
modeled for (A) 1980-1992 and (B) 2003-2009 periods, and 
differences between these two distributions (C). The agreement 
between presence-absence predictions (i.e., habitat suitability) was 
measured by summing the 90 predictions (threshold x iteration) for 
each reach of the Garonne river basin for each period, with color 
scale varying from green (no predicted presence) to red (90 predicted 
presences). The differences in the spatial distribution of the species 
were expressed with a color scale varying from blue (90 presences 
predicted only for 1980-1992) to red (90 presences predicted only for 
2003-2009). 



2704 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



I. Paz-Vinas ef al. 



Genetic and Demographic Analyses for Conservation 



■o 

CD 

'a. 



c 

CD 



30 



20 



10 



-10 



CD 

rr 



-20 



-30 — ' 



P< 0.01 



1980-1992 



1 

2003-2009 



Figure 5. Boxplots of the length of the occupied network by 
Parachondrostoma toxostoma in the Garonne river basin modeled for 
the periods 1980-1992 and 2003-2009. The length of the occupied 
network was the residuals of a linear regression linking the length of 
occupied network in both periods with the threshold setting method 
effect. 

likely (i.e., the first mill weirs date back from the 12th 
century, Blanchet et al. 2010). It is noteworthy that the 
bottlenecks highlighted here are "species-specific" rather 
than "basin-specific", given that for four other sympatric 
cyprinid fish species (i.e., Squalius cephalus, Leuciscus 
burdigalensis, Gobio gobio, and Phoxinus phoxinus), 
Paz-Vinas et al. (2013) demonstrated that bottlenecks 
were older (approximately 2000-6000 years ago) and of 
different magnitudes than those detected for P. toxostoma. 
We can hence reasonably conclude that the bottlenecks 
inferred here occurred during the last millennium and 
affected specifically P. toxostoma populations. 

Descriptive analyses revealed low levels of genetic diver- 
sity for all populations. Indeed, all diversity indices were 
up to approximately 3.3 times lower than those calculated 
for populations of other cyprinid fish species co-occurring 
with P. toxostoma in the Garonne river basin (Blanchet 
et al. 2010). They were all also remarkably lower than 
those calculated for P. toxostoma populations from the 



Rhone river basin (see Dubut et al. 2010). As an example, 
some microsatellite markers were monomorphic in cer- 
tain populations, whereas these same markers were highly 
polymorphic in populations from the Rhone river basin 
(Dubut et al. 2010). Similarly, Costedoat et al. (2005) 
demonstrated that the diversity measured at mitochon- 
drial genes for P. toxostoma was also significantly lower in 
the Garonne river basin than in the Rhone river basin, a 
result that may be a consequence of the recent coloniza- 
tion of the Garonne river basin from the Rhone river 
basin (i.e., approximately 57,000 years ago, Costedoat 
et al. 2005). Although the relatively poor genetic diversity 
found in the Garonne river basin probably has an impor- 
tant phylogeographical basis (Costedoat et al. 2005), it 
may reflect the more recent (200-700 years ago) and 
severe bottlenecks that we detected. 

Finally, our PGS also highlighted that P. toxostoma 
populations in the Garonne river basin were relatively 
homogeneous from a genetic standpoint. Indeed, most 
populations formed a single cluster with relatively low 
genetic differentiation within this cluster. This result sug- 
gests that these populations constitute a single panmictic 
unit at the basin level. There were, however, two notice- 
able exceptions to this general pattern; HER and SAV 
were genetically differentiated from all other populations. 
These two populations also demonstrated the lowest 
contemporary N e values, the lowest genetic diversities 
(i.e., H e , H 0j and AR), and the strongest bottlenecks. 
Altogether, this indicates that these populations may be 
discriminated from others (1) because gene flow between 
these populations and others are weak; and/or (2) because 
genetic drift and inbreeding were particularly high in 
these populations, causing divergence from other popula- 
tions in the Garonne river basin. 

To summarize, PGS provided a precise description of 
the current genetic state of P. toxostoma populations from 
the Garonne river basin. Overall, these results clearly indi- 
cate that long-term management should integrate the fact 
that the evolutionary potential of the species in this geo- 
graphic area may be weak. 

What did we learn from demographic data? 

Using time series abundance data at twelve locations, we 
found an overall demographic decrease of P. toxostoma 
populations that occurred in the last three decades. Evi- 
dence of a demographic decrease was further supported by 
comparing the P. toxostoma occurrence at the basin scale 
between two periods (1980-1992 and 2003-2009). This 
analysis revealed a significant decrease in the distribution 
range of P. toxostoma, representing 13.1% of the 1980- 
1992's distribution. These results confirm that over the 
range of the species, there is a decreasing trend in abun- 
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dance (Crivelli 2006; Poulet et al. 2011). This decrease 
contrasts with the increase in occurrence, abundance, and 
density of several sympatric species at the French scale 
such as Barbus barbus or Gobio gobio (Daufresne and Boet 
2007; Poulet et al. 2011). Despite this range-wide trend, 
we showed that not all local populations were subjected to 
a significant demographic decrease, as some of them 
display no particular trends, and one population even 
showed a significant demographic increase. There was no 
clear spatial pattern regarding these site-specific trends 
(see Fig. 2D). However, such site-specific analysis provides 
a basis for further analyses exploring the regional and/or 
local causes of demographic trends in the Garonne river 
basin. Indeed, a comparison implying healthy versus non- 
healthy (from a demographic point of view) populations 
may highlight the leading environmental factors affecting 
the demography of this species. 

To summarize, DMPs provided insights into the demo- 
graphic dynamics and changes in the spatial distribution 
of P. toxostoma in the Garonne river basin, which indi- 
cates that this species is ecologically weakened in this 
area, and thus restoration plans should be engaged to 
ensure the persistence of populations. 

Synthesis, implications, and conclusions: The 
conservation gain of combining genetic and 
demographic data 

Synthesis 

The history of P. toxostoma in the Garonne river basin is 
relatively recent and began approximately 57,000 years 
ago, when it colonized the Garonne from the Rhone river 
basin (Costedoat et al. 2005). Our results suggest that 
populations exhibited relatively large long-term N e 
(approximately 5000-8000 individuals per population) 
until severe and recent (approximately 800 to 200 years 
ago) demographic collapses entailed N e of less than a few 
hundred (sometimes less than a dozen) individuals. This 
means that very small numbers of effective breeders are 
currently sustaining populations in the Garonne river 
basin. This history led to genetically impoverished P. tox- 
ostoma populations in the Garonne river basin. Although 
most populations are genetically homogeneous, these 
demographic collapses also led to local differentiation in 
the Garonne river basin. In a more recent timeframe (i.e., 
the last two decades), we showed that this species experi- 
enced a global decrease in census size (N c ) over the entire 
Garonne river basin, although that some populations 
remained demographically stable or even increased locally. 
This recent decrease in N c was accompanied by a signifi- 
cant reduction of its spatial distribution over the Garonne 
river basin. Because both N e and N c are reduced in these 



populations, P. toxostoma in the Garonne river basin is 
confronted with a combination of ecological and evolu- 
tionary extinction risks, which reinforces its status of 
vulnerable species in the IUCN red list, and supports the 
implementation of conservation plans. 

Implications 

Our results illustrate how combining genetic and demo- 
graphic approaches is useful to target and to prioritize 
conservation and management plans for endangered pop- 
ulations. A main weakness of our study resides in the few 
number of sampling points common to both temporal 
trend and genetic analyses. However, this fact may well be 
the standard for most studies focusing on rare and threa- 
tened species. We therefore provide recommendations 
considering two cases. In the first case, both demographic 
and genetic are available at the sampling site level. In this 
case, combining genetic and demographic approaches 
allows identifying priority populations as those (1) having 
the lowest genetic diversity and N e ; and (2) being sub- 
jected to a significant and recent decrease in N c . For 
instance, we identified the Hers River as a priority popu- 
lation as both genetic and demographic indices are weak. 
In this case, we propose conservation strategies involving 
a program of stocking from broodstock stemming from 
healthy populations, combined with the restoration of 
habitat and connectivity with other rivers. Healthy popu- 
lations are those with stable N c and higher N e (such as 
the Petite Barguelonne and Louge rivers). In the second 
case, only one of the two metrics is available at the sam- 
pling site level. In this case, prioritizing conservation 
plans is less straightforward. For instance, some popula- 
tions (e.g., the Vendinelle River) were subjected to a 
sharp decrease in N c in recent years, however, no data are 
yet available regarding genetic diversity and N e dynamics. 
In this case, managers can conduct a genetic monitoring 
of these populations to help clarify the populations' sta- 
tus. On the other hand, some populations (e.g., the SAV) 
have low N e and low genetic diversity, but lack temporal 
data regarding N c . In this case, it is impossible to get the 
temporal trend of the populations. Thus, invoking the 
precautionary principle, we propose considering these 
populations as conservation priority. 

Conclusion 

To conclude, we showed how combining analyses based 
on point genetic studies (PGSs) and DMPs (i.e., a 
"demo-genetic approach") reveal complementary infor- 
mation underlying different processes operating at differ- 
ent timescales. Demo-genetic approaches allow (1) 
identification of "at risk" populations; (2) prioritizing 
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conservation and management actions; and (3) proposing 
plans that account for the evolutionary history and poten- 
tial of populations. We hence argue that demo-genetic 
approaches should be the norm in conservation practices. 
Indeed, these surveys would allow not only prioritizing 
and initiation of conservation plans (this study), but 
would also allow the evaluation of dispersal and connec- 
tivity through the use of genetic-based inference methods 
(Broquet and Petit 2009), as well as evaluation of the 
effectiveness of conservation plans (Schwartz et al. 2007; 
Osborne et al. 2012). We hope that this study will moti- 
vate conservation ecologists to invest in genetic monitor- 
ing, and conversely, conservation geneticists to initiate 
long-term demographic surveys. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Appendix SI. Current spatial distribution and recent 
distribution changes. 

Figure SI. Map of the Garonne river basin (South-Western 
France) representing (1) sites where P. toxostoma was 
unsuccessfully sampled for genetic analyses (white circles); 
and (2) sites where P. toxostoma was successfully sampled 
for genetic analyses (green circles). 

Figure S2. Maps representing sites where the occurrence of 
P. toxostoma was recorded (A) from 1 to 19 times during 
the 1980-1992 period and (B) from 1 to 14 times during 



the 2003-2009 period. 

Table SI. Information on microsatellite loci and multi- 
plexed PCR used in this study. 

Table S2. Observed number of alleles (na), expected (H e ) 
and observed (H 0 ) heterozygosities and departures from 
Hardy-Weinberg equilibrium (Vis) for all loci and popula- 
tions of P. toxostoma. No significant departures from 
Hardy-Weinberg equilibrium were found after applying 
Benjamini and Hochberg (1995) false discovery rate 
corrections. 

Table S3. Mean observed (H D ) and expected (H e ) hetero- 
zygosities, mean number of alleles over loci (NA), allelic 
richness (AR 8 for a minimum sample size of 8 individuals; 
AR 18 for a minimum sample size of 18 individuals) and 
departures from Hardy-Weinberg equilibrium (P is ) for each 
P. toxostoma population. 

Table S4. P-values for the linkage disequilibrium test for 
each pair of loci and population. 

Table S5. Results for the Wilcoxon's sign rank test com- 
puted by BOTTLENECK for each river and for the TPM 
and SMM microsatellite mutation models. 
Table S6. Median, 5 and 95% quartile values calculated for 
N 0 (the current effective population size), Ni (the past 
effective population size), T a (the time of the beginning of 
the demographic change, in years backwards from the pres- 
ent) and Log 10 (Nq/Ni; the magnitude of the demographic 
change) for each river, through the posterior distributions 
obtained with MSVAR 1.3. 

Table S7. Values for the Mann-Kendall's S statistic, vari- 
ance in S (Var[S]), mean densities and P-values obtained 
for the twelve time series with the modified Mann-Kendall 
trend test. 
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